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PREDICTING  ADSORPTION  EQUILIBRIA  OF  MIXTURES: 
COMPARISON  OF  POTENTIAL  THEORY  AND 
IDEAL  ADSORPTION  SOLUTION  THEORY 


1.  INTRODUCTION 

The  Dubinin-Radushkevich  (DR)  Equation  is  useful  for  correlating,  as  well  as 
predicting  pure-component  adsorption  equilibrium  data  for  adsorbates  that  are  liquids  at  the 
adsorption  temperature.  It  is  useful  because  (1)  it  is  a  relatively  simple  equation  that  fits 
adsorption  equilibria  data  of  non-polar  organic  adsorbates  reasonably  well  over  a  practical  range 
of  concentrations  and  temperatures,  (2)  the  theory  on  which  it  is  based  provides  a  means  of 
predicting  adsorption  equilibrium  using  liquid  properties  of  the  adsorbate,1  and  (3)  the  theory  can 
easily  be  extended  to  predict  the  adsorption  equilibrium  of  mixtures.2,3  A  second  approach  used 
to  predict  adsorption  equilibria  of  adsorbed  mixtures  from  pure-component  data  is  the  Ideal 
Adsorbed  Solution  Theory  (LAST)  of  Myers  and  Prausnitz.4  The  LAST  is  rarely  used  with  the 
DR  Equation.5  This  work  will  show  that  the  DR  Equation  can  be  successfully  incorporated  into 
the  LAST  framework. 

The  DR  Equation  is  based  on  a  theory  developed  by  Polanyi6  where  he  proposes 
that  the  important  intermolecular  interactions  controlling  the  adsorption  process  are  electronic 
dispersion  forces.  Further,  he  states  that  the  magnitude  of  the  electronic  dispersion  forces  is  both 
independent  of  temperature  and  proportional  to  the  molecular  polarizability  of  the  adsorbate. 

The  adsorption  potential  is,  likewise,  independent  of  temperature.  In  addition  to  the  discussion 
of  intermolecular  interactions  and  the  adsorption  potential,  Polanyi  notes  that  in  the  limit  of  low 
surface  coverage,  adsorption  isotherms  should  be  linear.  In  more  recent  literature,  this  is 
commonly  referred  to  as  Henry’s  Law.1,4,7,8,9,10 

Building  on  Polanyi's  work,  Dubinin11  found  that  many  adsorbates  have  potential 
curves  of  similar  shape,  and  he  proposed  a  relationship  having  two  adjustable  parameters  to  fit 
these  data.  This  is  the  basis  of  the  well-known  Dubinin-Radushkevich  (DR)  Equation.  Polanyi 
also  proposes  approximating  the  molar  polarizability  using  the  molar  volume,  thus  providing  a 
means  of  predicting  the  parameters  of  the  DR  Equation  from  readily  available  physical  properties 
of  the  adsorbate.  The  DR  Equation  has  proved  to  be  very  useful  for  making  practical  estimates 
of  adsorption  equilibria,1,12,13'14’15,16  but  it  is  often  criticized  for  lacking  a  Henry’s  Law  limit.10,17 

The  practical  advantage  of  the  Polanyi-Dubinin  theory  approach  is  that  it  is 
predictive.  If  adsorption  equilibria  data  for  an  adsorbate-adsorbent  pair  is  available,  then  the 
adsorption  isotherm  for  a  similar  adsorbent  on  the  same  adsorbent  can  be  estimated  knowing  the 
ratio  of  the  molar  volumes  of  the  two  adsorbates.18  In  keeping  with  Dubinin's  original  work, 
estimates  have  also  been  made  using  a  ratio  of  a  function  of  molar  polarizability.19  A  second 
advantage  is  that  the  same  theory  can  be  used  to  predict  adsorption  equilibria  of  mixtures.  The 
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theory  has  some  limitations.  Perhaps  the  most  significant  is  that  the  DR  Equation  and  its 
relatives  do  not  possess  a  Henry’s  Law  limit. 

Based  in  part  on  work  reported  by  Lewis,  et  al.,2  Grant  and  Manes3  discuss  a 
correlation  for  predicting  equilibria  of  mixed  vapor  adsorption.  This  approach  is  based  on 
Potential  Theory  principles  only.  These  principles  are  discussed  in  the  theory  section  of  this 
paper.  The  paper  of  Myers  and  Prausnitz4  also  outlines  a  theory  useful  for  predicting  adsorption 
of  mixtures.  The  IAST  is  well-tested  and  often  makes  good  predictions.1,4,9’1415  Regardless  of  its 
widespread  use  and  acceptance,  the  IAST  is  not  without  problems.  As  noted  by  Sircar  and 
Myers,7  the  reference  state  (the  spreading  pressure)  can  become  undefined  when  applied  with 
some  systems.  This  problem  is  especially  pronounced  near  saturation.7 

This  work  compares  the  multicomponent  Potential  Theory  approach  of  Grant  and 
Manes3  with  the  IAST.  Using  the  DR  Equation  to  correlate  pure-component  equilibria 
adsorption,  the  IAST  is  compared  mathematically  to  the  results  obtained  from  multicomponent 
Potential  Theory.  It  is  shown  that  the  IAST  and  Potential  Theory  treatments  of  multicomponent 
adsorption  produce  identical  results  when  the  correlation  coefficient  of  the  DR  Equation  is 
proportional  to  the  molar  volume  of  the  adsorbate.  This  is  an  assumption  of  the  original  DR 
Equation  but  has  been  discarded  by  some  researchers.  The  significance  of  the  Gibbs  adsorption 
isotherm  and  the  relationship  between  the  Gibbs  adsorption  isotherm  and  the  basic  assumptions 
of  the  DR  Equation  are  discussed.  It  is  shown  that  unless  the  pure-component  adsorption 
equilibria  relationships  satisfy  the  Gibbs  adsorption  isotherm,  the  IAST  cannot  be  applied 
because  the  resulting  equations  have  no  solution.  Further,  when  the  Gibbs  adsorption  isotherm  is 
satisfied,  it  is  shown  that  the  multicomponent  Potential  Theory  and  the  IAST  are  equivalent. 

2.  ideal  Adsorbed  solution  theory 

The  development  of  the  IAST  follows  closely  Gibbs’20  1906  treatment  of 
adsorption  on  the  surface  of  a  liquid.  Myers  and  Prausnitz4  point  out  that  thermodynamic  theory 
sets  limits  on  the  number  of  degrees  of  freedom  of  a  system.  The  phase  rule  for  adsorption 
proposed  is 

(degrees  of  freedom)  =  (number  of  components)  -  (number  of  phases)  +  3  (1) 

where  the  adsorbent  is  considered  "thermodynamically  inert"  (that  is,  it  is  not  considered  a 
component  itself).  The  phase  rule  for  adsorption  specifies  an  additional  degree  of  freedom  over 
the  phase  rule  for  bulk  solid-liquid-vapor  systems.  This  degree  of  freedom  implies  the  existence 

of  an  independent  thermodynamic  function  of  state.  In  IAST,  this  is  the  spreading  pressure,  it. 
Myers  and  Prausnitz  write 

Pi  =  *,•  pf(it) 


(2) 


where  p°  is  the  pure-component  adsorbate  partial  pressure,  xt  is  the  mole  fraction  of  the  ith 
component  in  the  adsorbed  phase,  and  ;ris  the  spreading  pressure  of  the  mixture.  The  spreading 
pressure  is  defined  by  the  Gibbs  adsorption  isotherm 


where  A  is  the  adsorption  area.  It  is  assumed  that  no  change  in  adsorption  area  will  occur  with 
mixing,  so  one  may  express  A  in  terms  of  pure-component  areas  and  adsorbed  phase  mole 
fractions: 


A=A1°xi  +  A2° x2  +  A3°x3...  (constant  T and  it)  (4) 

Since  the  adsorbed-phase  area  is  inversely  proportional  to  the  adsorbed-phase  concentration,  one 
may  write 


(5) 


where  n,  is  the  total  number  of  moles  adsorbed.  The  n°  are  pure-component  adsorbed-phase 
concentrations  which  correspond  to  the  spreading  pressure  of  the  mixture.  The  n°  are  expressed 
as  implicit  functions  of  the  spreading  pressure  by  Equation  3.  As  is  well  documented,  the  IAST 
makes  good  predictions,  even  in  some  cases  where  the  adsorbates  would  not  be  expected  to  form 
ideal  solutions  in  the  liquid  phase.4 

An  example  of  the  IAST  using  Langmuir  pure  component  isotherms  is  now  given 
to  emphasize  a  point  that  was  not  emphasized  by  Myers  and  Prausnitz;  the  theory  places  a 
restriction  on  the  pure-components  isotherms.  The  single  component  Langmuir  isotherms  are 

0 

In  p,  =  In  t — ! — r  -  lnfc,  (6a) 

1  (i-«i) 

In  j>2  =  In  .  - 1 n*;  (6b) 

V-Oi) 

where  fractional  surface  coverage  =  n  J  mx  is  the  adsorbed  concentration,  and  m  is  the  mono- 
layer  capacity.  Equations  6a  and  6b  are  inserted  into  Equation  3  to  determine  the  relationship 
between  the  adsorbed  phase  concentration  and  the  spreading  pressure  of  the  mixture. 
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nA 

RT 


(7) 


Equation  7  can  be  valid  only  if  0°  is  equal  to  02  and  mx  is  equal  to  m2.  For  pure- 
component  Langmuir  isotherms,  equal  mono-layer  coverages  are  required  for  thermodynamic 
consistency9  and  are  consistent  with  the  site-specific  assumption  used  in  deriving  the  Langmuir 
equation.  Using  the  equal  mono-layer  coverage  restriction  demonstrated  in  Equation  7  (all  n° 
are  equal)  one  may  apply  the  IAST  simply  by  multiplying  Equation  5  by  n*n,°  to  give  the 
following: 


ni=x\nt  +x2nt  =  «,  +n2  (g) 

Another  important  conclusion  that  can  be  drawn  from  Equation  7  is  that  &  is  a  unique  function 
Of  7T. 

3.  POTENTIAL  THEORY 

The  Potential  Theoiy  is  in  part  a  thermodynamic  theory;  but,  in  addition  to 
thermodynamics,  it  considers  the  nature  of  the  intermolecular  forces  acting  in  the  adsorption 
process.6  The  adsorption  potential  is  defined  as  the  difference  between  the  chemical  potential 
(the  chemical  potential  was  defined  by  Gibbs20)  of  a  substance  in  the  adsorbed  state  and  its 
chemical  potential  in  a  state  free  from  the  influence  of  an  adsorbent.  The  state  the  adsorbed 
substance  assumes  uninfluenced  by  the  adsorbent  is  taken  as  the  reference  state: 

e  =  n-firef  (9) 

In  addition  to  the  thermodynamics,  it  is  assumed  that  the  adsorption  potential  can 
be  expressed  by  an  equation  in  the  following  form: 

s  =  pF(0)  (10) 

The  adsorbent  contribution  F(0)  varies  depending  on  the  fraction  of  the  adsorption  space  filled 
0,  while  the  adsorbate  contribution  is  characterized  by  a  constant  /?.  Equation  10  follows  from 
the  assumption  regarding  the  nature  of  intermolecular  interactions  important  in  the  adsorption 
process:  the  interactions  are  due  to  electronic  dispersion. 

Building  on  the  work  of  Lewis,  et  al.,2  Grant  and  Manes  apply  the  Potential 
Theory  to  treat  multicomponent  adsorption.  Because  F{0)  of  Equation  10  is  independent  of 
properties  of  the  adsorbate,  the  ratio  sip  must  also  be  independent,  and,  therefore,  the  same  for 
all  adsorbates  of  the  mixture.3  Grant  and  Manes,  much  like  Myers  and  Prausnitz,4  list  a  series  of 
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equations,  which  when  solved,  simultaneously  (usually  numerically)  predict  multicomponent 
equilibrium  from  unconstrained  and  not  necessarily  related  pure-component  isotherms.  Grant 
and  Manes  briefly  discuss  the  conditions  under  which  their  theory  is  equivalent  to  the  LAST  but 
do  not  mention  any  required  relationship  between  pure-component  isotherms  used  with  their 
theory.  Belfort21  discusses  similarities  between  the  Grant  and  Manes  and  LAST,  but  due  to  an 
error  in  a  derivation,  reached  an  incorrect  conclusion  that  the  theories  were  equivalent  regardless 
of  relationships  between  pure-component  isotherm  equations. 

4.  MULTICOMPONENT  (MIXTURE)  EQUILIBRIA  ADSORPTION  USING  THE 

LAST  WITH  PURE-COMPONENT  DR  EQUATIONS 

Even  though  the  DR  Equation  lacks  a  Henry’s  Law  limit,  it  is  still  possible  to  use 
it  with  the  IAST.17  As  an  example,  the  DR  Equation  for  pure-component  adsorption  equilibria  is 
applied  to  the  IAST  pure-component  relationships.  The  DR  Equation  can  be  expressed  in  terms 
of  adsorbate  partial  pressure: 


In/?, 

+  In  p? 

(Ha) 

ln/?2 

+  In  pf 

(lib) 

The  fraction  of  the  adsorbed  space  filled  is  related  to  the  spreading  pressure  by  substituting 
Equations  11a  and  lib  into  Equation  3  to  get 


%A 

Jt 


■J-ln0,° 


f0d  J-lnff  = 


M 

RT 


V-ta02° 

fedyj-  In# 


00 


(12) 


where  the  fraction  of  the  adsorption  space  filled,  Ot,  =  nj  n-at,  n  is  the  adsorbed  concentration, 
and  nsat  is  the  adsorbent  capacity.  The  subscript  refers  to  the  component  number.  As  with  the 
Langmuir  equation,  the  function  of  6  in  the  integral  is  non-linear,  and  in  order  to  be  valid,  all  6° 
must  be  equal  (6?  =02°).  This  is  consistent  with  the  assumptions  of  the  Potential  Theory  for 
mixtures  and  similar  to  the  restriction  imposed  on  the  Langmuir  pure  component  isotherms  in 
Equation  7.  Equation  12  now  requires  the  following  relationship  be  true. 


A  =  *£ 
A 


(13) 


The  adsorbed-phase  concentration  at  saturation  is  defined  as 
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(14) 


where  W0  is  a  constant  representing  the  capacity  of  a  given  adsorbent  and  is  independent  of  the 
adsorbate  and  Vm  is  the  molar  volume  of  the  adsorbed  phase.  Equation  14  requires  that  the  ratio 
of  the  correlative  coefficients,  /?,  and  f32,  are  equal  to  the  ratio  of  the  adsorbed  phase  molar 
volumes  VmI  and  Vm2.  Unlike  the  case  with  the  Langmuir  isotherm,  the  saturation  capacity  of 
adsorbates  can  differ  for  each  pure-component  DR  Equation.  This,  too,  is  consistent  with  the 
theory  on  which  the  equations  are  based.  In  place  of  a  site-specific  assumption,  the  Potential 
Theory  assumes  the  adsorbate  fills  an  adsorption  space  and  the  amount  of  space  filled  is  related 
to  the  molar  volume  of  the  adsorbed  phase. 

The  conclusions  drawn  from  Equation  12  make  completing  the  IAST  derivation  using  the 
pure  component  DR  Equations  straightforward.  Multiplying  both  sides  of  Equation  5  by  n, 
gives  the  following  result: 


1  = 


(15) 


Dividing  denominator  and  numerator  of  each  term  of  the  sum  of  Equation  15  by  nsat,  results  in 


+ 

e?  el 


(16) 


and  since  0,°  =  one  may  derive  the  following  relationship: 

0°=  0X  +  02  (17) 

Thus,  as  expected  for  an  ideal  solution,  the  fraction  of  the  adsorption  space  filled  is  equal  to  the 
sum  of  the  fractional  filling  of  the  components  of  the  mixture. 

The  general  form  of  the  multi-component  DR  derived  from  DR  pure-component 

isotherms  is 


a- *,/>;•« p 


V.,E 

~RT 


(18) 


where  the  summation  is  over  ally  components,  including  component  i.  This  relation  gives  the 
adsorbate  partial  pressure  in  terms  of  the  adsorbed-phase  concentration  and  the  temperature. 

The  multi-component  form  of  the  DR  Equation  given  by  Yang1,22  is 


12 


(19) 


where  n,  is  the  total  adsorbed-phase  concentration,  W0  is  the  adsorbent  capacity,  Vmi  is  the 
adsorbed-phase  volume  of  the  rth  component,  k  =  \IE2  is  the  constant  that  characterizes  the 
adsorbent,  fa,  is  the  correlative  coefficient  of  the  rth  adsorbate  and  the  mole  fractions  in  the 
adsorbed  phase  are  related  by  Sx, •  =  1.  Algebraic  manipulation  reveals  that  Equation  19  and 
Equation  18  are  identical  if  $  =  Vm  i.  Since  Equation  13  shows  that  indeed  fat  =  Vmi,  then  one  can 
conclude  that  the  multi-component  formulation  of  Yang  and  the  LAST  are  equivalent  as  long  as 

0,°  =  02°. 


An  important  result  obtained  from  both  the  Langmuir  and  the  DR  derivations  is 
that  it  is  possible  to  determine  p°  as  a  function  of  &  of  the  mixture  as  opposed  to  n .  We, 
therefore,  can  write  Equation  2  in  the  equivalent  form: 


Pi 


(20) 


5.  THERMODYNAMIC  CONSISTENCY 


>otential  Theory  and  the  Gibbs  Adsorption  Isotherm. 


The  Gibbs  adsorption  isotherm  written  for  pure-component  adsorption  (Equation 
5)  plays  an  important  role  in  the  LAST.  This  equation  requires  that  there  be  a  relationship 
between  pure-component  isotherms  in  order  for  them  to  be  used  successfully  in  making  LAST 
predictions.  The  relationship  between  the  LAST  and  the  Potential  Theory  can  be  illustrated 
further  by  expressing  it  in  Potential  Theory  notation  and  evaluating  the  integral.  Making  the 
appropriate  substitutions,  Equation  5  becomes 

£p-p° 

[xA}=W0  ]0d{ft )  (21) 

Vo 

where  0  =  fa ini  /  W°  and  sx  =  -  RT In  — .  It  has  been  shown  that  for  non-linear 

Xj  p 

adsorption  equilibria,  the  spreading  pressure  must  be  a  unique  function  of  the  fraction  of  the 
adsorption  space  filled.  This  is  the  case  when  the  differential,  d[ Jf-) ,  of  Equation  21  is  of  the 

form  d  F{0).  Integrating  Equation  21,  keeping  in  mind  that  T is  held  constant,  results  in  the 
following: 
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s,  =  /3,F{8\T)  +  p.G(T) 


(22) 


If  the  potential  is  assumed  to  be  independent  of  temperature  and  if  G  =  0, 
Equation  22  is  identical  to  Equation  10  of  the  Potential  Theory;  thus  the  Potential  Theory  is 
consistent  with  the  Gibbs  adsorption  isotherm.  Furthermore,  when  using  the  IAST,  the  functions 
F{Q,T)  and  G{T)  must  be  the  same  functions  for  all  pure  components  of  the  mixture.  If  they  are 
not,  the  pure-component  isotherms  are  mutually  inconsistent  with  the  Gibbs  adsorption  isotherm, 
and  there  will  be  a  region  of  the  concentration  space  where  the  spreading  pressure  will  become 
undefined. 

5.2  Henry’s  Law  for  the  DR  Equation. 

The  DR  Equation  is  criticized  for  lacking  a  Henry’s  Law  limit.17  Heruy’s  Law  for 
adsorption  requires  that  in  the  limit  as  the  vapor-phase  concentration  approaches  zero,  adsorption 
isotherms  must  be  linear.  The  Langmuir  equation  has  a  Henry’s  Law  limit.  In  this  section  the 
limiting  behavior  of  the  DR  Equation  for  ideal  mixtures  is  explored.  Parameters  for  a 
hypothetical  binary  DR  Equation  are  shown  in  the  Table  below.  Several  Henry’s  Law  plots  of 
predictions  of  this  binary  equation  are  shown  in  the  Figure.  Whenever  some  of  Component  2  is 
adsorbed,  the  predicted  isotherm  of  Component  1  becomes  linear  as  its  concentration  approaches 
zero.  Thus,  the  Henry's  Law  is  satisfied  at  all  times  as  long  as  at  least  two  components  are 
adsorbed. 


Table.  Mixture  Parameters  For  The  DR,  Vm  —  J3. 


Cmipwieitt 

W  (npcg) 

1 

500 

100 

100 

2 

500 

200 

100 

6.  DISCUSSION 

It  is  important  to  note  that  when  F  is  independent  of  temperature  and  G  =  0,  the 
basis  of  the  Polanyi  Potential  Theory,  Equation  22  is  identical  to  Equation  10.  The  Potential 
Theory,  therefore,  conforms  to  the  restrictions  imposed  by  the  Gibbs  adsorption  isotherm.  The 
Gibbs  adsorption  isotherm,  however,  is  less  restrictive  than  the  Potential  Theory  and  allows  also 
for  a  function  of  temperature.  This  is  because  the  IAST  does  not  assume  anything  regarding  the 
nature  of  the  intermolecular  interactions  involved  in  adsorption. 

The  pure-component  DR  Equation  does  not  have  a  Henry’s  Law  limits  but  the 
multicomponent  DR  does.  If  there  is  more  than  one  adsorbate,  and  Component  2  is  maintained 
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Figure.  Natural  Henry's  Law  Limit  of  the  IAST-DR  Equation  for  a  Binary 
System  Consisting  of  Components  1  and  2. 

at  a  finite  concentration  while  the  concentration  of  Component  1  is  reduced  to  an  infinitesimal 
level,  Component  1 -Component  1  interactions  will  become  insignificant.  Component  2- 
Component  1  interactions  give  rise  to  a  Henry’s  Law  limit,  much  as  Component  1-adsorbent 
interactions  give  rise  to  a  Henry’s  Law  limit  for  the  Langmuir. 

The  Gibbs  adsorption  can  be  used  to  check  the  mixture  data  for  consistency. 
When  doing  LAST  calculations,  Equation  3  is  used  to  determine  the  pure-component  adsorbed- 
phase  concentration  corresponding  to  the  spreading  pressure  of  the  mixture.  When  there  is  a 
unique  spreading  pressure  of  the  mixture,  the  theory  places  a  condition  on  pure-component 
isotherms.  Therefore,  in  addition  to  mixture  isotherms,  two  sets  of  pure-component  isotherm 
data  can  be  checked  for  consistency  with  the  Gibbs  adsorption  isotherm. 

An  important  practical  consequence  of  this  is  when  Equation  22  is  satisfied  by 
analytical  pure-component  isotherm  equations,  the  IAST-multicomponent  form,  Equation  22,  is 
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also  analytical.  In  many  cases  when  pure-component  isotherms  are  used  with  the  IAST, 
numerical  solutions  are  unnecessary. 

7.  CONCLUSION 

The  IAST  and  the  Potential  Theory  mixing  rules  are  reviewed  with  particular 
emphasis  on  application  with  the  DR  Equation.  It  is  shown  that  the  IAST  and  the  Potential 
Theoiy  mixing  rules  either  differ  or  are  the  same,  depending  on  the  relationship  between  pure- 
component  equations  with  which  they  are  applied.  When  Equation  22  is  satisfied,  the  two 
approaches  give  identical  results.  When  Equation  22  is  not  satisfied,  results  differ.  Furthermore, 
when  Equation  22  is  not  satisfied,  the  numerical  method  does  not  converge. 

When  the  DR  Equation  is  used  with  the  LAST  to  derive  an  equation  to  describe 
mixture  adsorption,  the  resulting  equation  has  a  Henry’s  Law  limit,  as  well  as  a  finite  enthalpy  of 
adsorption  at  low  loading  of  a  given  component  i  as  long  as  another  component  is  present. 

When  mutually  consistent  pure-component  isotherm  equations  can  be  solved 
analytically  for  partial  pressure  in  terms  of  adsorbed-phase  concentration  and  temperature, 
mixture  equations  derived  using  the  IAST  are  also  analytical. 
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GLOSSARY 


E  adsorbent  potential  field  parameter  of  DR  equation  ($E  has  units  of  J/mol) 

F  general  function 

G  general  function 

p  partial  pressure  (Pa) 

p°  pure-component  partial  pressure  (Pa) 

psat  partial  pressure  at  liquid  saturation  pressure  (Pa) 

n  adsorbed-phase  concentration  (mol/kg) 

n°  adsorbed-phase  concentration  corresponding  to  the  pure-component  partial  pressure 
(mol/kg) 

nsa‘  adsorbed-phase  concentration  in  equilibrium  with  the  liquid  saturation  pressure  (mol/kg) 

n,  total  adsorbed-phase  concentration  of  a  mixture  (mol/kg) 

R  gas  constant  (J/mol-K) 

T  temperature  (K)  » 

Vm  adsorbed-phase  molar  volume  (cm3/mol) 

W0  adsorbent  capacity  {WJ$  has  units  of  mol/kg) 

x  mole  fraction 

y  variable 

Greek 

p  correlating  coefficient  of  DR  equation  (f \E  has  units  of  J/mol) 

s  the  adsorption  potential  (J/mol) 

tzA  product  of  spreading  pressure  and  adsorption  area  (J/mol) 

0  fraction  of  the  adsorption  space  filled 

0°  fraction  of  the  adsorption  space  filled  corresponding  to  the  pure-component  partial 

pressure 

Subscript 

i  property  of  the  ith  component 
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